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ABSTRACT 



We investigate stability of Toomre-Hayashi model for a self-gravitating, rotating 
gas disk. The rotation velocity, v v , and sound speed, c s , are spatially constant in the 
model. We show that the model is unstable against an axisymmetric perturbation 
irrespectively of the values of v v and c s . When the ratio, v^/c s , is large, the model disk 
is geometrically thin and unstable against an axisymmetric perturbation having a short 



radial wavelength, i.e., unstable against fragmentation. When v^/cs is smaller than 
1.20, it is unstable against total collapse. Toomre-Hayashi model of 1.7 < v v /c s < 3 
was thought to be stable against axisymmetric perturbations in earlier studies in which 



only radial motion was taken into account. Thus, the instability of Toomre-Hayashi 
model having a medium Vp/cg is due to not change in the surface density but that in 
the disk height. We also find that the singular isothermal perturbation is stable against 
non-spherical peturbations. 

Subject headings: gravitation — hydrodynamics — instabilities — stars: formation 



1. INTRODUCTION 

Toomre (1982) and Hayashi, Narita, & Miyama (1982) found an analytic model for a self- 
gravitating gas cloud. The model describes a disk in which the self-gravity balances with the 
centrifugal force and gas pressure. The rotation velocity, v v , and the sound speed, c s , are assumed 
to be constant. The density distribution of the model is expressed as 
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where 
and 

7 = 1 + |, (3) 

in the spherical polar coordinates, (r, 9, ip). 

This simple model is useful for studying dynamics of self-gravitating gas disk. Earlier re- 
searchers (see, e.g., Kalnajs 1971; Lemos, Kalnajs, & Lynden-Bell 1991; Syer & Tremaine 1996; 
Evans & Read 1998a,b; Goodman & Evans 1999; Shu et al. 2000) studied further simplified 
model, i.e., the razor-thin power law disk model, which is equivalent to the vertically averaged 
Toomre-Hayashi model. Besides of its simplicity, it is a good approximation to a centrifugally sup- 
ported disk formed by collapse of a rotating gas cloud (see, e.g., Matsumoto, Hanawa, & Nakamura 
1997; Saigo & Hanawa 1998). The centrifugally supported disk will evolve into single or binary 
stars associated with gas disks. The evolution may be driven by instabilities of the disk. Since the 
Toomre-Hayashi model is an exact equilibrium solution, it is ideal for linear stability analysis; we 
need not worry about inaccuracy of the unperturbed state. 

In this paper we concentrate on the simplest mode, the axisymmetric one. We will show that 
the Toomre-Hayashi model is unstable for any 7, i.e., any v^/cg. This result is different from 
prediction by Hayashi et al. (1982) and earlier analyses by Lemos et al. (1991) and Shu et 
al. (2000). While two-dimensional perturbation is taken into account in our analysis, only radial 
perturbation was considered in earlier analysis. We believe that our analysis will provide insights 
for stability of self-gravitating disk. 

The rest of our paper has the following structure. In §2 we derive energy principle for axisym- 
metric perturbation around the Toomre-Hayashi disk. In §3 we obtain the marginally stable mode 
for the Toomre-Hayashi disk. In §4 we show that it is unstable for any 7 and discuss implications 
of our stability analysis. Main conclusions are summarized in §5. 



2. PERTURBATION EQUATION AND ENERGY PRINCIPLE 

We consider a small perturbation around the Toomre-Hayashi disk in the polar coordinates. 
In this paper we use the displacement vector, £ = (£ r , £#), in the ip = const, plane to denote a 
perturbation. From the equation of continuity we can evaluate the change in the density as 

\ d 1 d 

Pi = - — 7T ( r2 P^r) — a ( sin #Me) , (4) 

where po denotes the density in the equilibrium and is given by Equation (1). From the angular 
momentum conservation we can evaluate the change in the rotation velocity as 
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where v^q denotes the rotation velocity in equilibrium. The perturbation equation for the equation 
of motion is expressed as 



and 
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8t 2 r tan 9 r 89 \po J ' r 89 
for the r- and ^-directions, respectively. The change in the gravitational potential, 0i, is related to 
that in the density by 
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Before detailed manipulation we consider the characteristic of the Tooomre-Hayashi model. 
Since the Toomre-Hayashi model has neither inner nor outer boundary, its global analysis is not 
straight forward (see Goodman & Evans 1999, for mathematical discussions on the stability of 
scale- free disks). 

The Tooomre-Hayashi disk is scale free and has no characteristic length scale. It has neither 
characteristic timescale; both the rotation period and sound crossing time are proportional to the 
radial distance, r. Since the disk has no characteristic radius, inner and outer parts of a disk are 
equally either stable or unstable. This means that the disk has no discrete eigenmode which grows 
exponentially in time (see also Goodman & Evans 1999). 

If the Toomre-Hayashi disk is unstable, the perturbation grows faster in an inner part and 
propagate outwards since the dynamical timescale is shorter at a smaller radius. To show this we 
derive the energy principle, 

T + W + S = (9) 
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r = 



w 



s = 



Po 



a JO 
b 



' 8t 2 



0% 



8t 2 



r 2 sin 9 drd9 , 



2M 2 + c s 



ro /*7T I 

Ja JO [ 

/*7T 

/ r 2 s'm9C (c 2 pi + p 4>i) d9 
Jo 



+ p\4>\\ r 2 sin 9 drd9 , 



(10) 

(11) 
(12) 



Equation (9) is derived from equations (4) through (8) by multiplying the complex conjugate of p\, 
£ r , and <f>i. We can set the inner and outer radii, a and b, arbitrarily. Equation (9) denotes the 
energy conservation in the second order in the interval of a < r < b. 

The term, T, denotes the product of acceleration and displacement. If it is positive, the 
perturbation grows and the system is unstable. The term, W, denotes the excess of rotation 
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energy, thermal energy and gravitational energy. We consider the condition that the term, S, 
vanishes. This condition is equivalent to 



&(r, 9) 
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Pl (r,6) = Pl (9)r ik - 5 / 2 
Mr, 9) = M0)r tk - 1/2 



(13) 
(14) 
(15) 
(16) 
(17) 



where k denotes the wavenumber in the logarithmic scale (= d/dlnr). When k ln(6/a) = 2n, the 
third integral vanishes in equation (9). We call this the periodic boundary condition in the following. 
Our periodic boundary condition is equivalent to the surface density perturbation assumed by 
Kalnajs (1971) and others. 

When we apply the periodic boundary condition, the term S is proportional to In (b/a) and 
scale free. Also the term T should be scale free and the growth timescale should be proportional to 
the radius, d/dt oc 1/r. This means that there exists no exponentially growing mode satisfying the 
periodic boundary condition. We seek, instead, a marginally stable mode satisfying the periodic 
boundary condition. 

Substituting the boundary condition into the perturbation equation we obtain 
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v v0 = C s ^2(7-1) . 

We solve equations (18) through (22) as an eigenvalue problem for k. 
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3. MARGINALLY STABLE MODE 

3.1. Singular Isothermal Sphere 

When 7 = 1, the Toomre-Hayashi disk is spherically symmetric and reduces to the singular 
isothermal sphere. The density distribution is expressed as 

* - si? ■ (25> 

Then the perturbation should be expressed by the spherical harmonics, Yp{9, ip). We concentrate 
on the axisymmetric perturbation for simplicity. Then the density perturbation is expressed as 

Pl (6) = AP e (cos9) , (26) 

where A and Pe denotes an arbitrary scale factor and the Legendre's polynomial, respectively. Then 
the change in the potential is given by 



4irGp 



1 1 - 1 

k 2 + -+£(£ + 1) 



(27) 



Substituting equation (27) into equation (20) we obtain 



k 2 = 7 --£(£ + l). (28) 

From equation (28) we conclude that the marginally stable mode exists only for £ = and the 
critical wavenumber is k = y/7/2. 

The critical wavenumber coincides with that of the equilibrium isothermal sphere. As noted 
in the textbook of Chandrasekhar (1939), the density distribution is approximately given by 
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(29) 



in the region ofr > c s / \J4irG Pc , where p c denotes the central density. This equilibrium isothermal 
sphere has almost the same density distribution as that of the singular isothermal sphere in the 
region ofr > c s /^4irGp c . The small difference between them can be regarded as a neutral 
perturbation on the singular isothermal sphere. Note the wavenumber, y/7/2, in Equation (29). 

The singular isothermal sphere is stable against non-spherical perturbation of £ ^ 0. Equation 
(28) leaves formally a possibility that the singular isothermal disk is unstable for any k. This 
possibility is, however, very unlikely. When k is very large, the wavelength is so short to destabilizes 
the gas sphere. Thus we can exclude the possibility that the singular isothermal sphere is unstable. 
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3.2. Toomre-Hayashi Disk 



In this subsection we solve the perturbation equation to obtain the critical wavenumber k as 
a function of 7. First we obtain 



89 



+ tan I 



(ik - v vl 



0, 



from equations (19) and (20). Equation (30) has the solution, 

Vtpl oc (sin 9) ik ~ 1 / 2 _ 
To avoid the divergence of v^i at 9 = 0, we seek modes of u^i = 0. 
Substituting v^i = into equations (18) through (22), we obtain 
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We solve Equation (35) with the boundary conditions d(f>i/d9 = at 9 = and ir as an 
eigenvalue problem. Once <\>\ is given, other perturbations, £ r , ^0, and pi are derived from Equations 
(32) through (34). For a given set of 7 and fc, we integrated Equation (35) with the Runge-Kutta 
method. By try and error, we found pairs of 7 and k for which Equation (35) has a solution 
satisfying the boundary condition at 9 = and ir. 

Figure 1 denotes the eigenvalue, k, as a function of 7. The solid curve connects the crosses 
which denotes the numerical obtained eigen values. The dashed curve does k = 7 for comparison. 
It is a good approximation to the numerically obtained eigenvalues for 7 > 1. 

Figure 2 shows P0/P1 as a function of 9 for the marginally stable modes. Each curve is labeled 
with 7. The relative density perturbation, pi/po, is symmetric with respect to 9 = tt/2 and has no 
zero point. When 7 is very large, the density perturbation is restricted near 9 = tt/2. Note that 
Figure 2 shows not p\ but pi/ po- The density perturbation, pi, is very small near the pole (9 = 
and 7r) even for 7 = 2, since po is very small there. 

After having completing numerical computations, we found an asymptotic solution for Equa- 
tion (35). To derive the asymptotic solution, we rewrite Equation (35) into 
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When 7 3> 1, </>i has an appreciable amplitude only near the disk plane (i.e., 8 ~ tt/2 or equiva- 
lently ( — 0). Thus, Equation (36) can be approximated as 



Equation (37) has an analytic solution 
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for k 2 = 7 2 — 1/4. This denotes an asymptotic solution in the limit of 7 » 1. 



Figures 3 and 4 show the marginally stable perturbation by cross section. Figure 3 shows the 
density distribution in the Toomre-Hayashi model of 7 = 10. The solid curves are the contours of 
log p = n/4 where n is an integer. Only the sector near the disk plane is shown in Figure 3. Figure 
4 is the same as Figure 3 but the marginally stable mode is superimposed on the Toomre-Hayashi 
model. As shown in Figure 4 the neutrally stable mode deforms the Toomre-Hayashi model into 
rings. 



4. DISCUSSION 

4.1. Instability Strip 

We have shown in the previous section that the Toomre-Hayashi model has a marginally 
stable mode for any 7. This means that the Toomre-Hayashi model of any 7 is unstable against 
axisymmetric perturbation, since a marginally stable mode is adjacent to an unstable mode. 

As shown by Hayashi et al. (1982), the Toomre-Hayashi model is unstable for total collapse 
when 7 < 2.441. In other words, the Toomre-Hayashi model of 7 = 2.441 has a marginally stable 
mode of k = 0. This marginally stable point is also shown by an asterisk in Figure 1. It is also 
shown in the previous section that the Toomre-Hayashi model of 7 = (i.e., the singular isothermal 
sphere) is unstable when k < y/7 /2. Accordingly, we conclude that the Toomre-Hayashi model is 
unstable in the lower right side of the marginally stable curve in Figure 1. 

The above conclusion apparently contradicts with the stability analyses by Lemos et al. (1991) 
and Shu et al. (2000) and prospect by Hayashi et al. (1982). All of them concluded that the 
Toomre-Hayashi model of 2.5 < 7 < 6 is stable against axisymmetric perturbation. 

The dash-dotted curves of Figure 1 denote the points of marginally stability obtained by 
Shu et al. (2000). They analyzed the stability of the Toomre-Hayashi disk under the thin disk 
approximation while taking account of magnetic fields. The magnetic field reduces the effective 
gravity and increases the sound speed (see, e.g., Li & Shu 1997; Shu & Li 1997; Nakamura & 
Hanawa 1997). The effects can be fully taken into account by replacing the gravitational constant 
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(G) and sound speed (c s ) with the effective ones. Thus the stability of the magnetized disk is 
essentially the same as that of non-magnetized one. The condition for the marginal stability of the 
axisymmetric mode is denoted by 

(g + 1/4) [Kik) - 1] 
7 2 - (P + 1/4) K{k) ' 1 ' 

where 

K(k) I r[(i/2 + ifc)/2]r[(i/2-ifc)/2] 

l j 2 r [(3/2 + ik)/2] r [(3/2 - ik)/2] ' 1 j 

denotes the Kalnajs function for m = 0. Lemos et al. (1991) obtained also the same dispersion 
relation although Equation (39) was not explicitly printed. 

As shown in Figure 1, the disturbance is oscillatory and stable when 7 is medium. This 
discrepancy is due to incompleteness in their search for unstable modes. Shu et al. (2000) 
considered only radial perturbations neglecting vertical structure. However, vertical perturbation 
is dominant in the marginally stable mode that we found in this paper. Thus they could not find 
this marginally stable mode and adjacent unstable mode. 

Equation (39) has an asymptotic form of 

7 = fe + 2 + ^ + 0(fc- 2 ), (41) 

in the limit of k 3> 3. Note that the solid and dash-dotted curves are parallel in the region of large 
k and hence for large 7. Change in the vertical structure has a significant effect even for a disk of 
a large 7. Comparison of our analysis and that of Shu et al. (2000) tells us importance of vertical 
perturbation in the instability of the Toomre-Hayashi model. The importance will be also true for 
stability of self-gravitating disk in general. 

Hayashi et al. (1982) prospected on the basis of their numerical simulations that the Toomre- 
Hayashi is stable when 7 is medium. Their search was also incomplete since the the dynamic range 
in the radial direction is restricted. They could not find an unstable mode having a small k. Note 
that the wavenumber of the marginally stable mode is 3.102 and 5.058 for 7 = 3 and 5, respectively. 
This means that the radius increases by a factor of 7.58 and 3.46, respectively, at every node of 
the marginally stable mode. Thus large dynamic range is necessary to find an unstable mode from 
numerical simulations. We think that we need dynamic range covering at least two wavelengths. 

As shown in the previous section, the change in the rotation velocity (5v v ) vanishes in the 
marginally stable mode. This means that the displacement is only in the vertical direction. The 
marginally stable mode changes hight of the disk but not the surface density. 



4.2. Eigenmode 



Although we have tried, we could not find the lower limit on k for a given 7. In other words, 
we failed in finding the other side of the instability strip. This failure is not due to lack in our 
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numerical technic. 

One might think that the other marginally stable mode would stem from the edge of the total 
collapse, (k, 7) = (0, 0.2441). It, however, cannot stem. The marginally stable total collapse has 
radial dependence different from those assumed in Equations (13) through (16). The change in the 
density distribution is expressed as 

« - t <42> 

^ - 2C tanh ( 7 C) 
oc r" 2 , (44) 



Po (43) 



for the marginally stable total collapse. This total collapse was called the breathing mode in Lemos 
et al. (1991). 

The other marginally stable mode might have different radial dependence which have been not 
yet known. Usually we assume some temporal and spatial dependences when we search unstable 
modes. Often we assume that the perturbation grows exponentially. An eigenmode may not grow 
exponentially in time in the Toomre-Hayashi disk. Remember our argument based on the energy 
principle in §2 (see also Goodman & Evans 1999, for another discussion on the stability of scale-free 
disks). 

The expansion wave solution of Shu (1977) gives us a hint. The eigenmode may have a form 

of 

K (45) 
The density distribution is expressed as 

"=2^(i)' <46> 

in the expansion wave solution. The difference from the singular isothermal sphere can be regarded 
as a perturbation propageting outwards while growing. As pointed out by Hunter (1977), the 
singular isothermal sphere can be regarded as a similarity solution. Similarly the Toomre-Hayashi 
model can be regarded as a similarity solution. An unstable mode can grow in proportion to the 
power of time when the unperturbed state is a similarity solution (see, e.g., Hanawa & Nakayama 
1997). 



5. CONCLUSION 



We discussed the stability of the Toomre-Hayashi model against axisymmetric perturbations 
in this paper. Our approach was conservative since we focused on the marginally stable mode. We 
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did not obtain an unstable mode. Nevertheless we had some new findings. This is because we used 
least approximations and took account of the vertical structure exactly. 

Although our analysis is restricted in the axisymmetric mode, the result will be helpful for 
more general analysis on the stability of self-gravitating disk. For later use we summarize the main 
results in the following. 

1. The Toomre-Hayashi model is unstable against axisymmetric perturbation irrespectively of 
its model parameter, 7. 

2. The wavelength of the unstable perturbation is proportional to the radius. The ratio of the 
wavelength to the radius is larger for smaller 7. 

3. Change in the vertical height can induce instability even when the surface density does not 
change. 

4. It is prospected that an unstable perturbation may grow in proportion of time in the power 
of time. 

5. The Toomre-Hayashi model of 7 = 0, i.e., the singular isothermal sphere is proven to be 
stable against a non-spherical perturbation. 

The authors thank Tetsuro Konishi for telling us the asymptotic solution shown in §3. This 
research is financially supported in part by the Grant-in-Aid for Scientific Research on Priority 
Areas and that for Encouragement of Young Scientists of the Ministry of Education, Science, 
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Fig. 1. — The critical wavenumber, k, is shown as a function of 7. The solid curve denotes the 
wavenumber of the marginally stable mode. The dashed curve denotes k = 7 for comparison. The 
dash-dotted curves labeled 'radial mode' denotes the the marginally stable mode obtained by Shu 
et al. (2000) under the thin disk approximation. 
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Fig. 2. — The neutrally stable mode is shown for the Toomre-Hayashi model of 7 = 2, 3, 5, and 
10. The ordinate denotes the relative density perturbation, pi/po while the abscissa does 9. 
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Fig. 4. — The same as Fig. 3 but the marginally stable mode is superimposed on the equilibrium 
model. 



